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Abstract 

We point out that the electroweak hne-tuning problem in the supersymmetric Standard 
Models (SSMs) is mainly due to the high energy definition of the fine-tuning measure. We 
propose super-natural supersymmetry which has an order one high energy fine-tuning measure 
automatically. The key point is that all the mass parameters in the SSMs arise from a single 
supersymmetry breaking parameter. In this paper, we show that there is no supersymmetry 
electroweak fine-tuning problem explicitly in the Minimal SSM (MSSM) with no-scale super¬ 
gravity and Giudice-Masiero (GM) mechanism. We demonstrate that the Z-boson mass, the 
supersymmteric Higgs mixing parameter /i at the unification scale, and the sparticle spectrum 
can be given as functions of the universal gaugino mass Mi/ 2 - Because the light stau is the light¬ 
est supersymmetric particle (LSP) in the no-scale MSSM, to preserve R parity, we introduce a 
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non-thermally generated axino as the LSP dark matter candidate. We estimate the lifetime of 
the light stau by calculating its 2-body and 3-body decays to the LSP axino for several values 
of axion decay constant fa, and hnd that the light stau has a lifetime in [10“^, 100] s for 
an fa range [10®, 10^^] GeV. We show that our next to the LSP stau solutions are consistent 
with all the current experimental constraints, including the sparticle mass bounds, B-physics 
bounds, Higgs mass, cosmological bounds, and the bounds on long-lived charge particles at the 
LHC. 



1 Introduction 


It is well-known that supersymmetry (SUSY) solves the gauge hierarchy problem in the Stan¬ 
dard Model (SM) naturally. From the low energy physics point of view, in the supersymmetric 
SMs (SSMs) with i?-parity conservation, gauge coupling unification can be achieved, the Light¬ 
est Supersymmetric Particle (LSP) such as neutralino can be a dark matter (DM) candidate, 
and the electroweak (EW) gauge symmetry can be broken radiatively due to the large top 
quark Yukawa coupling, etc. While from the high energy physics point of view, gauge coupling 
unihcation strongly suggests Grand Unihed Theories (GUTs), and only the superstring theory 
may describe the real world. Therefore, supersymmetry is a bridge between the low energy phe¬ 
nomenology and high-energy fundamental physics, and then is the most promising new physics 
beyond the SM. 

A SM-like Higgs boson with mass around 125 GeV was discovered in July 2012 from 
the hrst run of the LHG [1, 2]. This is indeed a little bit heavy in the Minimal SSM (MSSM) 
since it requires the multi-TeV top squarks with small mixing or TeV-scale top squarks with 
large mixing [3]. But it is hne in the next to MSSM (NMSSM) since the SM-like Higgs boson 
mass can be lifted further via the extra tree-level contribution and pushing up effect from 
the Higgs mass matrix diagonalization [4, 5, 6]. In addition, we have strong constraints on 
the parameter space in the SSMs from the LHG SUSY searches. For instance, the gluino 
mass rrtg and hrst two-generation squark mass rriq should be heavier than about 1.7 TeV if 
they are roughly degenerate rriq ~ rrig, and the gluino mass is heavier than about 1.3 TeV for 
rriq 3> rrig [7, 8]. Therefore, the naturalness in the SSMs is challenged from both the Higgs boson 
mass and the LHG supersymmetry searches. As we know, there are two kinds of dehnitions for 
hne-tuning measures: the low energy dehnition [9, 10, 11] and high energy dehnition [12, 13], 
which will be explained in details later. We shall point out that the naturalness conditions 
from the low energy dehnition can still be satished in principle, but the naturalness condition 
from the high energy dehnition is indeed a big challenge. However, supersymmetry is the 
connection between the low and the high energy physics, so we should consider seriously the 
hne-tuning problem from the high energy dehnition. Recently, it was found [14, 15] that the high 
energy hne-tuning measure [12, 13] will automatically be at the order one 0{1) in the J^-SU{5) 
models [16, 17, 18, 19] with the No-Scale supergravity (SUGRA) boundary conditions [20] 
and the Giudice-Masiero (GM) mechanism [21]. Based on the above results, we propose the 
super-natural supersymmetry whose high energy hne-tuning measure is at the order one 0{1) 
naturally. The essential point for super-natural supersymmetry is that all the mass parameters 
in the SSMs arise from a single supersymmetry breaking parameter. 

In this paper, we show that there is no residual electroweak hne-tuning (EWFT) left in 


2 


the MSSM if we employ the No-Scale supergravity boundary conditions [20] and GM mecha¬ 
nism [21], From no-scale supergravity, we have the supersymmetry breaking (SSB) soft terms: 
Mi /2 7 ^ 0 and mo = Aq = -Bo = 0, where M 1 / 2 , ^ 0 , ^0 and Bq are respectively the universal 
gaugino mass, scalar mass, trilinear A term, and Higgs mass mixing parameter at the GUT 
scale. We study various aspects of the No-Scale MSSM. We show how the parameter c, which 
is the ratio between the Higgs bilinear mass parameter ^{GUT) at the GUT scale and gaugino 
mass Mi/ 2 , changes as a function of the Z-boson mass Mz- In the no-scale MSSM, the light 
stau (ti) turns out to be the LSP, which has a dominant component. To solve the cosmolog¬ 
ical dark matter problem and keep the R parity intact, we introduce a non-thermally generated 
axino (a) as the LSP dark matter candidate. In other words, the light stau is the next to LSP 
(NLSP), and can decay to the LSP axino. The extremely weak interactions between axino 
and the lighter stau can cause the long-lived staus. The long-lived charged massive particles 
(sometimes also called as GHAMPs) like the stau in our case can disturb the successful predic¬ 
tions of Big Bang Nucleosynthesis (BBN), and may violate the other cosmological constraints 
such as catalyst BBN (GBBN) and structure formation bounds. Thus, we will estimate the 
lifetime of the NLSP stau by calculating 2-body {tr —>■ rd) and 3-body (tr —)■ rya) decays. 
In addition, in the viable parameter spaces which satisfy the sparticle mass bounds, B-physics 
bounds, and ruh = 125 ± 2 GeV, we obtain M 1/2 > 1800 GeV and tan/? > 27. We hnd that 
the NLSP stau can have lifetime from lO”"^ s to 100 s for the axion decay constant /„ range 
[10®, 10^^] GeV, which is consistent with all the above mentioned cosmological bounds. If we 
consider /a = 7 x 10® we estimate lifetime from 10“^ s to lO”"^ s depending upon jegl (will 
be dehned latter in the text) values, corresponding to the decay lengths of around a hundred 
meter to around 30 kilometers. Moreover, our NLSP stau satishes the recently reported ATLAS 
Gollaboration bounds on the long-lived staus as well. We also comment on the possible way 
of detection of staus in such a scenario. Furthermore, we present three benchmark points with 
gaugino masses about 1866 GeV, 2725 GeV, and 4589 GeV as well as Higgs boson masses about 
123 GeV, 125 GeV, and 127 GeV, respectively. The particle spectra are generically heavy due 
to large gaugino masses. 

This paper is organized as follows. In Section 2, we briefly review the No-Scale supergravity. 
In Section 3, we discuss the supersymmetry electroweak hne-tuning problems and propose the 
super-natural supersymmetry. In Section 4, we present the numerical calculations of EWFT and 
display various aspect of our study by presenting several graphs. In Section 5, we study dark 
matter in the No-Scale MSSM. We outline the detailed scanning procedure, and the relevant 
experimental constraints that we have considered. Results of our scans are given as well. A 

^Various astrophysical and cosmological considerations and axion searches suggest the lower limit of axion 
decay constant as /q > 6 x 10® GeV [22]. 
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summary and conclusion are given in Section 6 . 


2 Brief Review of No-Scale Supergravity 

Let us first briefly review the basic idea of No-Scale Supergravity, which was proposed to 
address the cosmological flatness problem [20] . It fulhll three constraints: i) the vacuum energy 
vanishes automatically due to the appropriate Kahler potential; ii) There exist flat directions 
that leave the gravitino mass M 3/2 undetermined at the minimum of the scalar potential (this 
is why it is called as No-Scale model); hi) The quantity StrA^^ is zero at the minimum. The 
large one-loop corrections would force M 3/2 to be either identically zero or of the Planck scale 
in case of violation of the third condition. A minimal Kahler potential that satishes the hrst 
two conditions is [ 20 ] 



( 1 ) 


K 


where T is a modulus held and <I>i are matter helds. They parametrize the non-compact 
SU{N, 1)/SU{N) X U{1) coset space. The third condition can always be satished in principle 
and is model dependent [23]. From the Kahler potential in Eq. (1) one automatically obtains 
the No-Scale boundary conditions mo = Ao = Bq = 0, while M 1/2 can be non-zero and 
evolve naturally, as is in fact required for supersymmetry breaking. The high-energy boundary 
condition Bq = 0 effectively Exes tan at low energy. This means the entire supersymmetric 
particle (sparticle) spectrum is determined by M 1/2 and in a very good approximation the whole 
sparticle spectra are linearly rescaled in terms of M 1 / 2 . The result is a natural one-parameter 
model, with M 1/2 the single degree of freedom for mass parameter. 

3 The Supersymmetry Electroweak Fine-Tuning Prob¬ 
lem and Super-Natural Supersymmetry 

To consider the hne-tuning issue in the supersymmetric SMs, we need to dehne the hne-tuning 
measures hrst. There are two kinds of definitions: the low energy definition and high energy 
dehnition. 

For low energy dehnition, there are two similar dehnitions. Let us review them one by one. 
The hrst dehnition considers the Higgs boson mass [9, 10]. The SM-like Higgs particle h in the 
MSSM is a linear combination of and H^. Its potential can simply be reduced to 


V = ml\h\^+ ^\h\^ , 


( 2 ) 
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where is negative. Minimizing the Higgs potential, we obtain the physical SM-like Higgs 
boson mass rrih 


ml = -2ml ■ 


The hne-tuning measure can be dehned as [9] 


Aft 


ml 


For a moderately large tan(3 = (H^)/{H^), for example, tan/? > 2, we have 


( 3 ) 

( 4 ) 


ml ~ |/i|2 + Itree + Irad , (5) 

where fi is the supersymmetric bilinear mass between and and and are 

the tree-level and radiative contributions to the soft supersymmetry-breaking mass squared for 
Hu- Therefore, we obtain the following naturalness bounds for 5% fine-tuning [10] 

• The p term or effective p term is smaller than 400 GeV. 

• The square root Mi = -|- of the sum of the two stop mass squares is smaller 

than 1.2 TeV. Consequently, we can show that the light sbottom mass is lighter than mi^. 

• The gluino mass is lighter than 1.8 TeV. 


The second low energy dehnition, which is similar, considers the Z boson mass. With the 
one-loop effective potential contributions to the tree-level MSSM Higgs potential, we get the 
Z-bosom mass Mz 

Ml _ {m]j^ + E;^) - {mlj^ + T.I) tan^ (3 _ 3 
2 tan^ /3 — 1 ^ ’ 

where and are the contributions arising from the one-loop effective potential defined in 
Ref. [11] and tan/? = ^ is the ratio of the two Higgs Vacuum Expectation values (VEVs). All 
parameters in Eq. ( 6 ) are dehned at the electroweak scale Mew- To measure the EWFT, we 
dehne [ 11 ] 

Ch, = WhJ {tav? 13 - 1 ) 1 , Chu 
~ 1 ) 1 ) C'su = 

The hne-tuning measure is dehned by [11] 

Aew = max( 2 Gfc)/(M|) . ( 8 ) 


: I — mjj^ tan^ /?/(tan^ /? — 1 ) |, = | — |, 

— E“tafo/?/(tafo/? — 1)1 . (7) 
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Note that Aew only depends on the weak-scale parameters of the SSMs, and then is hxed by 
the particle spectra. Hence, it is independent of how the SUSY particle masses arise. 

For high energy dehnition in the GUTs with gravity mediated supersymmetry breaking, the 
typical quantitative measure Aeenz for hne-tuning is dehned by the maximum of the logarithmic 
derivative of Mz with respect to all fundamental parameters a* at the GUT scale [12, 13] 


Aeenz = Max{AfUT} ^ aGUT 


d\n{Mz) 

91n(af^'^) 


( 9 ) 


Because the SM-like Higgs boson mass in the MSSM is smaller than the Z boson mass at 
tree level and can be lifted by the top squarks radiatively, the discovery of Higgs boson with 
mass around 125 GeV [1, 2] indeed constrains the viable parameter spaces in the SSMs. In 
the MSSM, we might need the large trilinear soft term At or say large stop mixing if we want 
stop masses around 1 TeV [3]. In the NMSSM, the SM-like Higgs mass can also be lifted 
via the additional tree-level contribution and pushing up effect from the Higgs mass matrix 
diagonalization [4, 5, 6]. So the NMSSM looks more natural. In short, the 125 GeV Higgs boson 
mass will not induce big hne-tuning in the SSMs from the low energy dehnitions of hne-tuning 
measures. On the other hand, although the null results of the LHG Runl have raised the lower 
bounds on the masses of gluino, hrst/second generation squarks, and sleptons [7, 8, 24, 25, 26], 
they are still within the upper bounds from the 5% hne-tuning requirements via the low energy 
dehnitions given above. Therefore, the SSMs are still hne if we allow a few percent hne-tuning 
from the low energy dehnitions. The key problem is the high energy dehnition of hne-tuning 
measure. For example, we can have the benchmark points which have the low energy hne-tuning 
measure Aew around 20 while the high energy hne-tuning measure Aeenz around 1,500. For 
example, see the benchmark points 1 and 2 in Table 1 of Ref. [27]. 

Because the hne-tuning measures for high energy dehnition in the viable SSMs are very large 
at the order of 10^ (0(10^)), we would like to explore the supersymmetry breaking scenario 
whose hne-tuning measure for high energy dehnition is automatically at the order one (0(1)). 
In other word, the hne-tuning measure in Eq. (9) is exactly one in the dream case. Interestingly, 
the solution with Aeenz = 1 is indeed simple. If there is one and only one mass parameter 
M* in the SSMs, Mz is a trivial function of M*, and we have the following approximate scale 
relation 


M” = U (g) M: , (10) 

where /„ is a dimensionless parameter, and q denote the dimensionless coupling parameters, 
such as gauge and Yukawa couplings, as well as the ratio between p and M 1/2 in the MSSM, 
etc. 
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For the nearly constant /„ of Eq. (10), we have 


f ^ ’ 

and therefore we obtain 

d\n.{M'^) dM^ M” 5M^ 1 

ain(Mf) “ ~ 'MjlMj ~ TJ"" 

Consequently, the hne-tuning measure is an order one constant 


(91n(M") 

(91n(Mp) 


~ 0 ( 1 ) . 


( 11 ) 


( 12 ) 


(13) 


Therefore, there is no electroweak hne-tuning problem in such kind of SSMs. In the no-scale 
SU{5) model where the fi problem is solved via the GM mechanism [21], we obtain /i cx M 1/2 cc 
M 3 / 2 , suggesting mutual proportionality, where the gravitino mass can be situated around 30 
TeV to elude the gravitino problem. Crucially, hnding /i ~ M 1/2 approximately rescales the 
sparticle spectra per variation only in M 1 / 2 . Taking M* = M 1 / 2 , it has been shown that there 
is indeed no residual EWFT problem [14, 15]. 

Based on the above discussions, we propose the super-natural supersymmetry with Aeenz — 
1. The necessary conditions for the super-natural supersymmetry are 


• The Kahler potential and superpotential can be calculated in principle or at least inspired 
from a fundamental theory such as string theory with suitable compactihcations. In 
other words, one cannot add arbitrary high-dimensional terms in the Kahler potential 
and superpotential. 

• There is one and only one chiral superheld or modulus which breaks supersymmetry. And 
all the supersymmetry breaking soft terms are obtained from the above Kahler potential 
and superpotential. 

• All the other mass parameters, if there exist like the fi term in the MSSM, must arise 
from supersymmetry breaking. 


Therefore, all the supersymmetry breaking soft terms and mass parameters in the SSMs are 
linearly proportional to the gravitino mass. 

In this paper we will show that there is no residual EWFT problem in the MSSM with the 
No-Scale SUGRA boundary conditions and GM mechanism. In order to achieve our goal, our 
general strategy is as follows. With the gauge coupling unihcation, we shall determine the GUT 
scale Mqut- Then we obtain the ratio /i/Mi/ 2 , top quark Yukawa coupling, and all remaining 
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input parameters at Mgut- As a result, Mz is a trivial function of M 1 / 2 . Therefore, from the 
above discussions, there is no electroweak hne-tuning problem in such No-Scale MSSM. We 
shall confirm this in the following via numerical calculations of Aeenz for n = 1 in Eq. (13) and 
point out that it is the same for the traditional choice of n = 2 or any other positive integer 
n since all the mass parameters are linearly proportional to M 1/2 in the No-Scale MSSM with 
the GM mechanism. 

4 The Numerical Calculations of Aeenz in the No-Scale 
MSSM 

We use publicly available code SuSpect2.43 [28] for our calculations. We employ the no-scale 
SUGRA boundary conditions (mo = Aq = 0) and generate points randomly for the following 
parameter space 

300 GeV <Mi /2 < 2000 GeV , 

l<tan/?<60. (14) 

We use sign(/i) > 0, m* = 173.3 GeV [29] and m^ = 4.25 GeV. We implement Bq = 0 
requirement numerically with a width |i?o| ^ 1 GeV that is comparable to the electroweak 
radiative corrections. In Fig. 1, we plot ^{GUT) and fine-tuning measure respectively on the 
left and the right vertical axes while M 1/2 is plotted on the horizontal axis. Let us discuss our 
calculations and Fig. 1 in more details. We divide our calculations into two parts. 

• Part 1: We vary M 1/2 and tan/3 in the interval given above Eq. (14) with no-scale SUGRA 
boundary conditions and demand \Bo\ < 1 GeV. The thick blue line in Fig. 1 is consisted of 
points we generated. Here we see that niGUT) oc M 1/2 or niGUT) = cMi /2 where c is a 
proportionality constant. We obtain the value of parameter c by htting our points with a hrst 
degree polynomial curve with no constant, which is shown as a thin blue line. The value of 
parameter c turns out to be 1.128, i.e., ^{GUT) = 1.128Mi/2- Note that the value of parameter 
c may depend on the distribution of points. Hence the points should be distributed as evenly 
as possible. We have tried our best to have evenly distributed points. 

• Part 2: Now we calculate hne-tuning using Eq. (12). Remember that here the deriva¬ 
tive is approximated by For the justihcation of the approximation, the variation 

of Mi /2 and Mz should be very small compared to their original values. We select the points 
for FT-calculations from evenly generated points in Part 1. We vary tan/3, M 1/2 and Mz- We 
insist on two requirements when we make the variation of M 1 / 2 , tan/3 and Mz'- 




1. The following GUT parameters vary as small as possible: Mgut, three Yukawa couplings 
HiiGUT) and three gauge couplings gk{GUT)). 

2. n{GUT) generated in program is as close to the value Hgut = I.I 28 M 1/2 as possible. 

For the first requirement, we first analyze the variation of GUT parameters when tan/3, M 1/2 
and Mz are varied and make sure that either they do not vary or the variation is negligible. 
For the second requirement, we set the criterion manually that the absolute value of relative 
change of/i(GUr), i.e., |(p(G'UT)'-/i(GUr))/p(G'UT)| = |(/i(GUT)'-1.128Mi/2)/1.128Mi/2| 
(where g,{GUT)' is the value of g,{GUT) in any given run when we vary M 1 / 2 , tan/3 and Mz ) 
should be less than 0.00001. Thus, we vary M 1/2 and Mz and collected several hundred points 
which fulhll the above criterion. Maroon points in Fig. 1 represent our results of fine-tuning 
measure calculations. We see that when the thick blue line (our data) and the thin blue line 
(the fitted line) are somewhat apart (at M 1/2 around 1000 GeV and 1800 GeV respectively), 
hne-tuning is also relatively large (~ 1.195 and 1.911 respectively). When these two lines meet, 
we see that hne-tuning measure is about 1 {Aeenz ~ !)■ 

In addition, we study the dependence of parameter c on the Z-boson mass Mz and present 
our results in Fig. 2. Here, we perform random scans by assuming tuq = Hq = 0 and Eq. (14) for 
individual values of Mz from 82 GeV to 102 GeV. Fig. 2 shows various values of the parameter 
c as a function of Mz- We see that there is a couple of irregular points in the curve at Mz= 
83 GeV and 86 GeV, respectively. In rest of the curve, we see that at Mz= 82 GeV, the value 
of c is around 1.23 and as the values of Mz increases, the value of c decreases until Mz= 89 
GeV. After that c rises with Mz for a short while until Mz= 93 GeV and starts decreasing 
again as Mz increases to 102 GeV. In Fig. 3, we present graphs in M 1/2 — tan/3 plane. Here we 
show that the initial values of tan/3 increases, though steadily, with M 1/2 for hxed Mz, as well 
as with Mz implicitly. We display plots only for six values of Mz to show the trend, although 
we have generated data, as discussed above for individual values of Mz from 82 GeV to 102 
GeV. Similarly, in Fig. 4, we present results in M^-tan/3 plane, where we select two different 
values of M 1 / 2 , 800 GeV and 1600 GeV, and pick up the corresponding tan /3 values along with 
corresponding Mz values. This hgure clearly shows the dependence of initial values of tan /3 
on Mz values. The plateau that we see for Mz values 90 GeV to 93 GeV corresponds to the 
values of the parameter c from 0.6 to 1.15. 

Moreover, we show sparticle masses as functions of M 1/2 in Fig. 5. In this figure, the value 
of parameter c = /i(GUr)/Mi /2 varies only slightly in the range c = 1.128 ± 0.001. In the left 
panel we have a plot for color sparticles. The green, purple, orange and black points represent 
rrig, mi^, rriu^ and respectively. In the right panel, the dark green points represent 
m-±, blue points represent m^o, red points show and brown points depict mgi- Thus, all 
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the sparticle masses are indeed linearly proportional to the gaugino mass. We would like to 
emphasize that our calculations of Aeenz for n=l is equivalent to the traditional choice of 
n=2 since all the mass parameters are linearly proportional to M 1/2 in the No-Scale MSSM 
with the GM mechanism. In fact, we generically have 

M” = (15) 

Therefore, |cIln(M^)/cIln(M(), 2 )| ~ 1 is valid for any positive integer n. 



W|/,[GeV] 


Figure 1: ^{GUT) (left-vertical axis) and EWFT measure Aeenz (right-vertical axis) versus 
Mi/ 2 - The thick blue line represents data generated, and the thin blue line is the line htted to 
the data. Maroon points depict the values of Aeenz- 


Figure 2: Parameter c versus Mz plot. 
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Figure 3: Plots in tan/3-Mi/2 plane. 
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Figure 4: Plot in the tan /3-Mz plane. 
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Figure 5: Sparticle masses as functions of Mi/ 2 - See text for color coding. 
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5 Dark Matter in No-Scale MSSM 


There are compelling evidences available in favor of the existence of dark matter (DM) [30]. 
Experiments suggest that it consists of electrically neutral stable massive particles [31, 32], In 
the MSSM with i?-parity conservation, the lightest neutralino (y?) can be the LSP and then 
the dark matter candidate [33, 34] (For reviews, see [35, 36].). The observed DM relic density 
can be obtained if we assume that neutralino pairs annihilate through pseudo-scalar Higgs 
(H-resonance), Z-boson (Z-resonance) or even through Higgs boson (h-resonance) to the SM 
particles in the early Universe. Also, the LSP neutralino can coannihilate with the NLSP, such 
as the lighter stau (ri), lighter stop (U), chargino (x]*"), etc, to reduce their number density 
and satisfy the observed relic density bounds. In the CMSSM/mSUGRA, even with these 
mechanisms available to satisfy the relic density bounds, there is still a large parameter space 
where the relic density of neutralino is unusually large or the LSP is a charged sparticle like 
the lighter stau. In this scenario, to preserve the i?-parity, one can assume that Xi D Eire 
the NLSPs and can decay to the non-MSSM sparticles. For example, both the axino a, the 
fermionic partner of the axion and the gravitino, the fermionic partner of the graviton can be 
the LSP candidates. 

In the no-scale MSSM, the LSP turns out to be the light stau, since the neutralino mass 
increases faster in RGE flow with M 1/2 as compared to the lighter stau. In Ref. [37], it was 
pointed out that the problem of a stau LSP can be alleviated when the unihcation scale is 
raised sufficiently above the GUT scale, leaving the possibility open for a bino LSP as the dark 
matter candidate 

In our case, the gravitino should be very heavy as discussed above, so we are left with 
the d LSP scenario. The axino is a promising dark matter candidate beyond the MSSM [39, 
40, 41, 42, 43, 44, 45, 46]. It appears when we extend the MSSM with the Peccei-Quinn 
mechanism [47, 48] in order to solve the strong GP problem. The mass of axino depends 
on the models and SUSY breaking scheme [49, 50, 51, 52]. The extremely weak interactions 
between the axino and charged slepton allow a long-lived NLSP charged slepton such as the 
lighter stau (ri). These long-lived staus, which were also called charged massive particles 
(GHAMPs) in Refs. [53, 54, 55, 56, 57, 58], can be probed at the LHG and have been studied 
extensively [59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71]. 

However, the Long-lived staus can cause problem and violate various cosmological bounds 
such as the Big Bang Nucleosyntheis bounds [72], catalyst BBN bounds [73, 74, 75] and struc¬ 
ture formation bounds [76]. Thus, in such a scenario, one should have a mechanism through 
which the lifetime of the lighter stau via decay to the axino LSP is within reasonable time, 

®See [38] and references therein for phenomenological issues related to No-Scale models. 
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and then the above mentioned cosmological bonnds can be escaped. Snch a scenario has been 
stndied assnming stan as the NLSP and gravitino as the LSP in Refs. [77, 78, 79, 80]. However, 
the gravitino mass shonld be heavier compared to the gangino mass in no-scale snpergravity, 
and the parameter space satisfying all the constraints discnssed above have heavy gangino mass 
Mi /2 ^ 1800 GeV from the later stndy. Thns, we can only assnme a as the LSP. In the follow¬ 
ing, we will explore the parameter space systematically where the light stan (fi) is the NLSP 
and the axino (a) is the LSP with the no-scale SUGRA bonndary conditions nnder varions 
collider and cosmological constraints. 

5.1 Phenomenological Constraints 

Using the no-scale SUGRA bonndary conditions, we perform random scans over the parameters 
space given in Eq. (14). Also, we consider /i > 0, = 173.3 GeV, and = 2.83 

GeV. The initial data are generated with SnSpect2.43 [28]. In this stndy we nse the following 
parameter space: 


300 GeV <Mi /2 < 5000 GeV , 

1 < tan < 60 . (16) 

The data points collected all satisfy the reqnirements of the radiative electroweak symmetry 
breaking (REWSB), no tachyonic sfermions or psendo-scalar Higgs bosons. We then employ 
MircOmegas3.6.9.2 [81] to calcnlate the relic density and other constraints. After collecting the 
data, we reqnire |i?o| ^ 1 GeV and the following bonnds on sparticle masses from the LEP2 
experiment 


mfj >105 GeV 
m^± > 103 GeV 

Also, we nse Ref. [82] to implement the following B-physics constrai 

10"^ < BR(R, ^ /i+/i-) < 6.2 X 10"^ (2a) 
10“^ < BR(6 ^ sy) < 3.87 x 10“^ (2a) 

0.15 < “ < 2-41 (3a) 


0.8 X 

2.99 X 


(17) 

(18) 
(19) 


( 20 ) 

( 21 ) 

( 22 ) 
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In addition, we impose the following bounds from the LHC as well 


nih = 123 - 127 GeV [1, 2 ] 
nig > 1.7 TeV (for irig ~ niq) [7, 8 ] 
rrig > 1.3 TeV (for rrig -C rriq) [7, 8 ] 


(23) 

(24) 
(26) 
(26) 


As far as the muon anomalous magnetic moment is concerned, we require that the benchmark 
points be at least as consistent with the data as the Standard Model. 

5.2 Numerical Results 

Here we present results of our scans. In the left panel of Fig. 6 , we show plot in Mi/ 2 -tan /3 plane 
for Mi /2 ^ 3000 GeV. The gray points satisfy the constraints of the REWSB and no tachyonic 
sfermions or pseudo-scalar Higgs bosons, blue points satisfy the sparticle mass bounds, green 
points form a subset of blue points and satisfy B-physics bounds, and red points form a subset 
of green points and are consistent with rrih = 125 ± 2 GeV constraint. In this plot we see that 
in the left top corner around tan/3 > 30, we do not have points with successful REWSB. We 
note that the application of the constraint |i3o| ^ 1 GeV reduces the parameter space to a strip 
of points like Fig. 3. We also notice that a sharp cut in blue points around M 1/2 ~ 1 TeV 
due to the requirement niq > 1800 GeV. When we apply the Higgs mass bounds (red points), 
in our present data we need M 1/2 > 1800 GeV and tan/3 > 27. As we discussed earlier, in 
no-scale MSSM, the gravitino would be very heavy due to M 3/2 oc M 1 / 2 . This clearly shows 
why we have chosen axino LSP scenario instead of gravitino LSP scenario. For comparison, in 
the right panel, we also display the same graph with M 1/2 ^ 5000 GeV. We see in this hgure 
that we have more red points with slight increase in tan /3 range. For our further analyses, we 
will restrict ourselves with the data M 1/2 < 3000 GeV. As we will show in the benchmark point 
Table 1 that even with M 1/2 < 3000 GeV, we already have heavy spartcile spectra. 

For our DM studies we assume non-thermally produced (NTP) axino LSP from decays of 
stau NLSP. We assume that each stau NSLP decays to an axino. The relic density of such 
axinos satisfy the observed DM relic density, ~ ^dm- As we mentioned earlier, because 

of extremely weak interaction with axino LSP, stau NLSP can be long-lived. Since the long- 
lived charged particles like stau can disturb successful predictions of the BBN and violate the 
other cosmological bounds, it is essential to calculate the lifetime of stau NLSP. This requires 
the knowledge of the axino mass. Usually the axino mass is a model dependent parameter. 
Note that axino is not a particle in the MSSM, we calculate its mass in the following way. As 
we just stated, we assume the observed dark matter relic density = 0.11 [32] and use the 
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following relation to compute the axino mass 

= 'rriNLSP-^ - 7^ ■ (27) 

In order to have mj^LSP (in our case the light stau is the NLSP) and fl^LSph"^, we present 
the results in — VtjsiLSph^ plane in Fig. 7. The gray points represent successful REWSB 
and no tachyonic sfermions or pseudo-scalar Higgs bosons. The aqua points satisfy sparticle 
mass bounds, B-physics and Higgs mass bounds. The purple points are subset of aqua points 
and satisfy |Ho| ^ 1 GeV constraint. Black solid line represent observed dark matter relic 
density while black dashed line show the upper WMAP9 5a bound. We pick up those points 
which have flNisph"^ greater than the WMAP9 5a bounds. This gives us both the lighter stau 
mass and the value of relic density to be used in Eq. (27). We show the NLSP stau mass and 
corresponding LSP axino mass in Fig. 8. Here, we see that the axino LSP mass range [565, 638] 
GeV corresponds to the stau NLSP mass range [788, 920] GeV. The stau NLSP (in our case the 
right-handed stau tr) can decay into the axino via tr —)■ ra (2-body decay), tr —)■ ray (3-body 
decay) and tr —)■ rdqq (4-body decay). The lifetime of stau NSLP is governed by the 2-body 
decay, while the contribution of 3-body decay remains below about 3%. The 4-body decay 
channels like tr —)■ rdqq, tr —)■ rayy and tr —)■ Tdl^l~ have negligible impact on the total stau 
decay rate as they are suppressed by an additional factor of a compared to the 3-body decay. 
This is why we do not consider them in our present study. The calculations of the decay rates 
T^tr —)■ ra) and P(rii; —)■ ray) have been reported in [86, 87] considering SUSY hadronic or 
KSVZ axion models [88, 89]. In [86], these calculations were done using an effective theory where 
(s)quark were integrated out, and then may lead to logarithmic divergences. It was also shown 
that these logarithmic divergences can be regulated with the cutoff fa, and only the dominant 
contributions are kept. In addition, a factor ^ and a mass scale m are introduced to parametrize 
the uncertainty associated with this cutoff procedure. We will use the decay rates given in [87] 
in our work because of absence of such issues. In the calculation of P(r^ —)■ ra), the scalar 
one-loop self-energy function (Hq) and vertex function (Gq) [90] are evaluated by the LoopTools 
package [91]. For the calculations of F(fH —)■ ra) we use Eqs. (3.1-3.3) and for F(rj? —)■ ray) 
we use Eqs. (3.4-3.6) in Ref. [87]. In numerical calculations we take a = a^^ijnz) = 1/128, 
sin^ 6w = 1 — rriwl^'z — 0 .2221, IcqI = 1/3 and y= 1, where [cqI and y are the electric charge 
and Yukawa coupling of the new heavy KSVZ quarks [87]. The lifetime (r^J of the light stau as 
a function of its mass for four different values of /„ is shown in Fig. 9. The green, orange, 
blue and red points represent solutions for the axion decay constants fa = 10®, 10^®, 10^^, 
and 10^^ GeV, respectively. The solid black line represents r?^ =1 second as a reference line. 
As expected, we see that the lifetimes increase with /„ values. For fa = 10® GeV, 10^^ GeV, 
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and 10^^ GeV, the lifetime of light stau is less than abont 10“^ second, 1 second, 100 seconds, 
respectively. In Fig. 14 of [87], plots in plane are shown with /„ = 10^^, 3 x 10^^, 10^^ 

GeV, \eQ\= 1/3 and mbino = 1.1 mf^,1.01 The anthors of Ref. [87] have given a very 
comprehensive presentation of varions cosmological constraints on stan NLSP and axino LSP 
masses. In Fig. 14 of [87], they have considered = Qdm and implemented strnctnre 

formations bonnds, primordial nncleosynthesis bonnds (both the hadronic and electromagnetic 
bonnds). In order to compare onr resnits with Fig. 14 of Ref. [87], we have to keep in view 
both Figs. 8 and 9 of onr paper. First of all we note that onr solntions lie in two narrow 
windows ~ [565, 638] GeV and ~ [788, 920 GeV]. From Fig. 14 of [87], we see that 
the strnctnre formation bonnds exclnde axino mass rria ^ 100 GeV for fa ~ 10^^ GeV and 
axino mass aronnd 250 GeV for fa = 10^^ GeV. Moreover, the hadronic BBN constraints 
(both the conservative and severe) imposed by the denterinm (D) abnndance mostly remain 
important for axino mass less than 100 GeV bnt rise np to higher axino mass ranges ~ 500 
GeV for fa = 10^^ GeV. Similarly, the corresponding electromagnetic BBN constraint remains 
important for stan mass less than 400 GeV where axino mass is in the range 1 GeV to 300 GeV 
for fa valnes 10^^ and 10^^ GeV. In scenario where stan NLSP has the lifetime greater than 
10 ^ s, negatively charged stan can form bonnd states with the primordial nnclei and thereby 
catalyse the formation of ®Li and ®Be in the early Universe [74], which is why they are called 
catalyse BBN (GBBN). From Fig. 14 of Ref. [87], we see that onr stan NLSP having lifetime 
> 10^ s may be rnled ont. As far as the stan NSLPs with lifetime Tf^ < 100 s is concerned, as 
also indicated in Ref. [87], one can expect the mild constraints as compared to above mention 
bonnds. We have not fonnd the detailed stndy discnssing the cosmological constraints for 
Ri ~ 100 s. When we compare onr resnits with the above observations, we note that onr 
solntions seem hne althongh there is a slight difference in assnmption (we do not have hxed 
the ratio between bino mass and stan mass). Maybe a separate analysis is needed to stndy 
the stan NLSP in no-scale MSSM with lifetime greater than 100 s with npdated astrophysical 
constraints. Note that Fig. 9 shows that depends on the fa and \eQ\. From above 
discnssions, it is evident that the cosmological bounds indeed restrict the parameter spaces of 
the SSMs with stau NLSP and axino LSP. Moreover, one can probe and study the properties 
of such axinos at colliders. As mentioned earlier, since axino interactions are extremely weak, 
the direct detection of axinos seems hopeless. Likewise, their direct production at colliders 
is very strongly suppressed. Instead, one expects a large sample of the NLSPs from the pair 
productions or cascade decays of heavy sparticles, provided the NLSP belongs to the MSSM 
spectrum. In our case, the stau NLSPs can be produced in cascade decays of heavy sparticles 
or from directly pair produced. One can measure the mass of staus using time-of-flight data 
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from the muon chambers. In addition, it is essential to study the decays of staus to measure 
and axino mass. Our analysis given above clearly shows that the lifetime of stau NLSP can 
be very large (small) depending upon the parameter space under consideration and parameters 
such as fa and \eQ\. We can measure the lifetimes of these charged massive particles at collider 
according to their decay lengths (cr). These lifetime measurements at colliders can be divided 
into three scenarios depending upon the decay length (cr^J and the size of the detector L 
which is typically L < 0.1 km. In the hrst scenario -C 0{cm), with very short lifetime 
-C In the second scenario the lifetime is supposed to be 0{cm) < < L, while 

in the third scenario we have > L. In Fig. 10, we present the lifetime (in the left panel) 
and the decay length (in the right panel) of staus for our data. Here we have used the minimal 
value of fa with two different values of jegl to see what are the minimal lifetimes and decay 
lengths in our present data. In this hgure, the aqua points correspond to fa = 7 x 10® and 
I eg I = 1/3, and the purple points represent fa = 7x 10® and jegl = 1. In the left panel, we see 
that for aqua and purple points the lifetimes are respectively between 10“® s to 10“^ s and 
10“^ s to 10“® s. In right panel, the minimal decay length for purple points (|eg| = 1) is less 
than 0.1 km, while for aqua points (|eg| = 1/3) it is about 5 km. This hgure shows that with 
heavy stau masses, the decay length of stau can be within the detector L < 0.1 km. If such 
staus decay to axinos, it is possible that we may observe LSP axino scenario with |eg| = 1. 
Moreover, it is estimated in Refs. [92, 93] that for a relatively small lifetime 10“® to 10“® s, 
we can still observe a substantial number of in-hight decays of staus within the detector. It 
is expected that the staus with large lifetimes may loose their energies through ionization and 
slowed down. These staus can be stopped within the main detector [94, 95, 96], or an additional 
dedicated stopping-detector can be used [97, 98, 99]. In this case, it is possible to study the 
subsequent decays of trapped staus by recording their stopping points/times via those devices. 
Furthermore, there is a proposal to place a water tank around the detector [100] or use the 
surrounding rock as a detector for the study of trapped stau decays [101]. In this paper we do 
not study the 2-body and 3-body decays of stau NLSP for our data at colliders and leave it 
as a future project. Such studies can be found in Ref. [87]. The collider phenomenology with 
the directly produced long-lived staus within the MSSM and calculations for dominant cross 
section contributions to stau pair productions are reported in Ref. [102]. In fact, the collider 
searches for charged long-lived particles have been performed at the LEP [103, 104, 105, 106] 
, HERA [107], the Tevatron [108, 109, 110], and the LHC [111, 112]. Recently, the ATLAS 
Collaboration [113] has reported that the long-lived stau, in the SSMs with gauge-mediated 
supersymmetry breaking (GMSB), is excluded up to its masses between 440 and 385 GeV for 
tan /3 between 10 and 50, with a 290 GeV limit in the case where only the direct stau production 
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is considered. In the context of the simplihed Lepto-SUSY models, where sleptons are stable 
and have a mass of 300 GeV, sqnark and glnino masses are exclnded up to a mass of 1500 
and 1360 GeV, respectively. In addition, the stau mass below 500 (339) GeV is excluded for 
the direct+indirect (direct only) production by the GMS Gollaboration [114]. While the LHGb 
Gollaboration has looked for stau with mass 124 and 309 GeV without any success and has put 
limits on di-stau production cross section [59]. Because the GMSB scenario is assumed in the 
latest studies of GHAPMs at the LHG, these stau mass bounds are for a different kind of the 
SSMs. The implications of these results on the No-Scale MSSM would necessitate a detailed 
study. At the moment, just taking these stau mass limits, it seems to us that our stau NLSP 
solutions are heavy enough to evade the LHG bounds. 

We display three benchmark points in Table 1 with gaugino masses about 1866 GeV, 
2725 GeV, and 4589 GeV as well as Higgs boson masses about 123 GeV, 125 GeV, and 127 GeV, 
respectively. Point 1 represents the scenario where we have data satisfying all the current con¬ 
straints discussed previously, as well as have M 1/2 < 2000 GeV and Higgs mass around 123 
GeV. Point 2 depicts the scenario where we have M 1/2 between 2000 GeV to 3000 GeV and 
Higgs mass mass around 125 GeV. Point 3 is an example of points with M 1/2 ^ 3000 GeV 
and Higgs mass around 127 GeV. We can see that in all these Points the sparticle spectra are 
heavy. For instance, for Point 1, we see that the masses for the hrst two families of squarks 
are slightly above 3000 GeV, while the stops are around 3000 GeV and glnino mass is around 
3895 GeV. Similarly, the masses for the left-handed slepton of the hrst two families are heavier 
than 1000 GeV but the right-handed slepton masses are about 679 GeV, while 2 are 587 
GeV and 1184 GeV, respectively. Electroweakinos are also heavy and in the range of [800,2000] 
GeV. For Points 2 and 3, the particle spectra are even more heavy. The masses for the hrst 
two families of squarks are above 4700 GeV and 7600 GeV respectively for Points 2 and 3. 
Similarly the light stops are about 4000 GeV and 6400 GeV. The left-handed slepton masses 
for the hrst two families are respectively 1739 GeV and 2890 GeV for Points 2 and 3 while 
the right-handed slepton masses are 987 GeV and 1653 GeV. The NLSP staus are 841 GeV 
and 1374 GeV for Points 2 and 3, respectively. Such heavy spectra can not be probed at the 
14 TeV LHG, which will provide a strong motivation for 33 TeV and 100 TeV proton-proton 
colliders. In Ref. [115], it was shown that the squarks/gluinos of 2.5 TeV, 3 TeV and 6 TeV 
may be probed by the LHG14, High Luminosity (HL)LHG14 and High Energy (HE) LHG33, 
respectively. This clearly shows that our models have testable predictions. Moreover, in the 
future if we have collider facility with even higher energy, we will be able to probe over even 
larger values of sparticle masses. 
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Figure 6: The viable parameter spaces in M 1/2 — tan/? plane. The gray points have the 
successful REWSB and no tachyonic sfermions or pseudo-scalar Higgs bosons. The blue points 
further satisfy sparticle mass bounds, green points form a subset of blue points and satisfy 
B-physics constraints, and red points are subset of green points and satisfy Higgs mass bounds. 
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Figure 7: versus The gray points have the successful REWSB and no tachyonic 

sfermions or pseudo-scalar Higgs bosons. The aqua points further satisfy the sparticle mass 
bounds, as well as the B-physics and Higgs mass bounds. And the purple points are subset 
of aqua points and satisfy \Bq\ < 1 GeV. The horizontal black solid line represents r2h^=0.11, 
while dashed black line shows 17/2^=0.137. 
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Figure 9: The lifetimes of staus (xfj in seconds (s) as functions of stau masses for |eQ|=l/3 
and y=l. Here, the light-green, orange, blue and red points represent the axion decay constants 
fa = 10®, 10^®, 10^^, and 10^® GeV, respectively. The horizontal black solid line represents 
= 1 second. 
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Figure 10: In the left panel, the lifetimes of staus (r^J in seconds as functions of stau 
masses for y=l and fa = 7x 10® GeV. The aqua and purple points respectively represents 
I eg I = 1/3 and jegl = 1. In the right panel, with the same color coding, corresponding decay 
lengths (cTfj) in kilometers (km) as functions of 


6 Discussions and Conclusion 


We emphasized that the electroweak hne-tuning problem in the SSMs is mainly because of the 
high energy dehnition of hne-tuning measure. We proposed the super-natural supersymmetry 
with the order one high energy hne-tuning measure naturally since all the mass parameters 
in the SSMs arise from a single supersymmetry breaking parameter. With the numerical cal¬ 
culations of the hne-tuning measures, we showed explicitly that we indeed do not have the 
supersymmetry electroweak hne-tuning problem in the MSSM with no-scale supergravity and 
Giudice-Masiero (GM) mechanism. We studied various aspects in the No-Scale MSSM, and 
found that the Z-boson mass, the fi{GUT) parameter and the sparticle spectra can be given as 
functions of the universal gaugino mass parameter Mi/ 2 - Introducing a non-thermally generated 
axino LSP, we estimated the lifetime of stau NSLP by calculating its 2-body body and 3-body 
decays to the axino LSP for various values of axion decay constant fa- In particular, our stau 
NSLP solutions satisfy all the current experimental bounds, such as sparticle mass bounds, B- 
physics bounds, Higgs mass bounds, cosmological bounds, and the bounds on long-lived charge 
particles at the LHG. 
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Point 1 

Point 2 

Point 3 

Afl/2 

1865.849 

2725.144 

4588.77 

mo 

0.0 

0.0 

0.0 

Aq 

0.0 

0.0 

0.0 

tan (5 

28.129 

30.594 

33.822 

Bq{GVT) 

0.599 

0.251 

-0.716 

/i{GUT) 

1991 

2788 

4431 

SqCewsb) 

92 

118 

4292 

/i(EWSB) 

1970 

2735 

166 

rrih 

123 

125 

127 

mu 

2054 

2816 

4321 

rriA 

2054 

2815 

4321 

mH± 

2056 

2816 

4322 

m^o 

820, 1520 

1219, 2238 

2100, 3802 

m^o 

1978, 1985 

2747, 2753 

4310, 4317 

m~± 

Xl.2 

1520, 1985 

2238, 2753 

3802, 4317 

nig 

3895 

5546 

9030 

'^UL,R 

3505, 3358 

4968, 4746 

8035, 7604 

^ h .2 

2786, 3225 

3957, 4554 

6411, 7343 

^dL.a 

3506, 3356 

4968, 4717 

8035, 7604 

mr 

01,2 

3205, 3247 

4530, 4574 

7268, 7364 

777/£>12 

1199 

1737 

2889 

rrii^a 

1177 

1701 

2819 

^eL,R 

1202, 679 

1739, 987 

2890, 1653 

m^i 2 

587,1184 

841, 1705 

1374, 2822 


Table 1: All the masses are in this table are in units of GeV and fi > 0. All this points satisfy 
the sparticle mass, B-physics constraints described in Section 5.1. Points 1, 2 and 3 represent 
parameter space with ruh = 123 GeV, 125 GeV and 127 GeV respectively. 
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